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ABSTRACT 



We test the recently published all-sky reddening map of Schlegel, Finkbeiner 
& Davis (1998 [SFD]) using the extinction study of a region in the Taurus 
dark cloud complex by Arce & Goodman (1999 [AG]). In their study, AG use 
four different techniques to measure the amount and structure of the extinction 
toward Taurus, and all four techniques agree very well. Thus we believe that 
the AG results are a truthful representation of the extinction in the region and 
can be used to test the reliability of the SFD reddening map. The results of our 
test show that the SFD all-sky reddening map, which is based on data from 
COBE/DIRBE and IRAS/ISSA, overestimates the reddening by a factor of 1.3 
to 1.5 in regions of smooth extinction with Ay > 0.5 mag. In some regions of 
steep extinction gradients the SFD map underestimates the reddening value, 
probably due to its low spatial resolution. We expect that the astronomical 
community will be using the SFD reddening map extensively. We offer this 
Letter as a cautionary note about using the SFD map in regions of high 
extinction {Ay > 0.5 mag), as it might not be giving accurate reddening values 
there. 

Subject headings: dust, extinction — infrared: ISM: continuum 
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1. Introduction 

Extinction and reddening caused by interstellar dust effects the detected emission from 
most astronomical sources in the sky. In most galactic and extragalactic studies, especially 
when studying so called "standard candles," the effects of dust on the source's detected 
brightness and color need to be taken into account. Hence it is very desirable to know the 
extinction (and reddening) anywhere on the sky. 

Recently Schlegel, Finkbeiner and Davis (1998; hereafter SFD) published an all-sky 
reddening map, based on satellite observations of far-infrared emission (at 100 and 240 
/im) from dust. This new reddening map will surely be used by many researchers seeking 
to apply reddening "corrections" for their work, so it is very important to verify its 
accuracy and reliability. The SFD map uses data obtained by the DIRBE (Diffuse InfraRed 
Background Experiment) on board CORE combined with ISSA (IRAS Sky Survey Atlas) 
images. The COBE/DIRBE experiment had better control of absolute calibration than did 
IRAS, but a larger beam (0.7°; as compared with ~ 5' for IRAS). SFD use the DIRBE 
data to calibrate the IRAS/ISSA images, and, after sophisticated processing, they obtain 
a full sky map of the 100 /im emission from dust, with point sources and the zodiacal 
light foreground removed.^ Their resulting reddening map is more accurate and has better 
resolution (~ 6.1') than the previous existing all sky reddening map of Burstein and Heiles 
(1978, 1982). 

In order to test the extinction map derived from the SFD analysis, we compare it with 
Arce & Goodman's (1999; hereafter AG) recent extinction study of a region of the Taurus 
dark cloud complex. The AG study uses four different techniques to measure extinction 



■^See SFD for a full description of their elaborate foreground zodiacal light subtraction 
technique, point source extraction and overall data reduction. 
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along two constant-Right Ascension "cuts" of several degrees in length. The four techniques 
in AG utilize: 1) the color excess of background stars for which spectral types are known; 
2) the ISSA 60 and 100 /im images; 3) star counts; and 4) an optical {V and R) version of 
the average color excess method used by Lada et al. (1994). The study finds that all four 
methods give generally similar results, and concludes that all four techniques are reliable 
ways to measure extinction in regions where Ay ^ 4 mag. In this Letter we compare the 
extinction map derived from SFD to the AG extinction map which is based on ISSA 60 and 
100 /im images in the Taurus region. This comparison provides a test of the reliability of 
the SFD reddening map for regions of extinction {Ay) higher than 0.5 mag. 

2. Comparing the two results 

In this section we compare the extinction derived from the SFD reddening map 
(hereafter Avgp^) to the extinction AG obtained using the ISSA 60 and 100 /im maps 
(hereafter Avjgg^)- In order to be consistent, we must compare both methods with the same 
spatial resolution. As part of their image reduction procedure, SFD convolve the ISSA 
images with a FWHM=3.2' Gaussian, which results in ISSA images with 6.1' resolution. 
For the purpose of comparing the two extinction maps, we also convolve the AG ISSA 
images with a FWHM=3.2' Gaussian, so that both extinction determinations have the same 
resolution. 

When AG compare Ay^gg^ with the extinction measured with star counting techniques 
and the average color excess method (techniques 3 and 4 above), they do so by plotting 
extinction versus declination. The extinction values shown in Figures 3 and 5 of AG are 
actually traces of the extinction averaged over the ~10' span (in the R.A. direction) of 
their CCD fields. For consistency we consider the same kind of average, and we obtain the 
SFD reddening {E{B — V)sfd) for the two 10' wide constant-R.A. cuts of AG and average 
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E{B — V)sFD over Right Ascension. In order to convert from the color excess value given 
by the SFD maps to visual extinction, we use the equation Ay = RvE{B — V), where Ry, 
the ratio of total-to-selective extinction, is equal to 3.1 (SFD; Kenyon et al. 1994; Vrba & 
Rydgren 1985). 

In Figure 1 we plot extinction versus declination for the two constant R.A. cuts in the 
Taurus dark cloud complex region. In addition to Ay^pj^ and Ay^^^^, we show the extinction 
traces of the extinction obtained from the average color excess and star count methods of 
AG, for reference. We will concetrate our discussion on the ISSA (AG) and SFD methods. 
It can be seen that these two methods trace similar structure in both cuts. The extinction 
peaks associated with the IRAS core Tau Ml (Wood et al. 1994), the dark cloud L1506, 
the dark cloud B216-217, and the IRAS cores Tau B5 and Tau Bll (Wood et al. 1994) are 
easily detected. The 1 a error (not including systematics) of Ay^^^^ is 0.12 mag (AG). SFD 
quote an uncertainty of 16% in their reddening estimates. For the most part, in regions 
unassociated with a pronounced peak in extinction, Ay^^^ is about a factor of 1.3 to 1.5 
larger than Ay^^^^. It is also important to note that in regions where there is a peak in the 
extinction, the value of ^V/ssa i^ closer to that of Ay^^^ than in the rest of the trace, and 
in some of the peaks Ay^^^^ exceeds Ay^^^. 

There are two places where the two traces appear to show very different extinction 
structures. In both of these places Ay^^^^ shows a steep dip: one is near the peak associated 
with B216-217 in cut 2 (around declination 26.9°); and the other is near the end of the cut, 
around declination 28.3° (see Figure 1). These two regions coincide with the position of 
IRAS point sources, and will be discussed in more detail further on. 
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3. Comparison between the two methods 

SFD use data from the IRAS and DIRBE experiments to construct a full-sky map of 
the Galactic dust based on its far-infrared emission. The IRAS data are used as a source 
of 100 /im flux images, and the DIRBE data are used for absolute calibration and as a 
source of 240 //m data. SFD derive the dust color temperature using the ratio of 100 to 
240 /im emission from DIRBE0, and then use this dust color temperature to convert their 
ISSA-based maps of 100 /xm emission to maps of dust column density. As a result of this 
procedure, the SFD ISSA-based 100 /im emission map has a spatial resolution of 6.1', but 
their dust color temperature map has a resolution of ~ 1.4° (Schlegel 1998). 

To transform from the map proportional to dust column (hereafter Tsfd) to a 
reddening {E{B~V)) map, SFD use the correlation between the intrinsic B — V color of 
elliptical galaxies and the Mg2 line strength. The Mg2 line strength of an elliptical galaxy 
correlates well with its intrinsic B — V color so that the Mg2 line index can be used along 
with photometric measurements of the galaxy in order to obtain a reasonably accurate 
measurement of its reddening (Faber et al. 1989). SFD use a procedure where the residual 
of the B — V color versus the Mg2 line index for 389 elliptical galaxies is correlated with the 
estimated reddening from their maps, using a Spearman Rank method. The resultant fit is 
then used by SFD to convert from tsfd to reddening {E{B — V)sfd) at each pixel. We 
then convert reddening to extinction, using the relation Ay^pj^ = RvE{B — V)sfd, with 
Rv = 3.1 (SFD). 

In order to obtain their ISSA-based extinction, Avj^g^, AG begin by using ISSA 60 



^Note that deriving dust temperature from a long wavelength color ratio, such as 100/240 
fim is superior to using the 60/100 //m ratio, because the effects of both point sources and 
transient heating of small grains are minimized. 
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and 100 /iin images to obtain a dust color temperature map from the flux ratio at each 
pixel. This color temperature is then used to convert the observed 100 /im flux to 100 /im 
optical depth, rioo- As in SFD, the next conversion, from dust opacity to extinction, is tied 
to a separate technique of obtaining extinction. AG chose to use a method similar to that 
described in Wood et al. (1994), which itself is ultimately based on work by Jarrett et 
al. (1989). Jarrett et al. (1989) correlate 60 ^m optical depth (tqo) with extinction (Ay) 
obtained from star counts, and Wood et al. (1994) multiply the Jarrett et al. tqo values by 
100/60 to derive a conversion from rioo to Ay, assuming low optical depth. Thus, using 
Wood et al.'s technique, AG's conversion of rioo to Ay ultimately rests on Jarrett et al.'s 
correlation of 60 /im optical depth with star counts. As described in AG, this tqq — Ay 
correlation is very tight for Ay < 4 mag, so we ascribe very little of the uncertainty in 
^VissA to this conversion. 

4. The cause of the discrepancy 

Arce & Goodman (1999) obtain the extinction toward the Taurus region using four 
different techniques. All four give similar results in terms of the absolute value and 
overall structure of the extinction. Most importantly, AG find that Ay^^^^ agrees well 
with extinction values obtained by measuring the color excess of background stars for 
which spectral types are known, which is the the most direct and accurate way to measure 
reddening (see Figure 2 in AG). So, the question now is why does the SFD method give 
extinctions which differ from the four determinations of extinction made by AG? 

Absolute values. Although data reduction and zero point determination of the emission 
is very elaborate and accurate in the SFD method, the normalization of the reddening per 
unit flux density (conversion between tsfd and E{B — V)sfd) seems to overestimate the 
reddening in regions of high dust opacity. As explained above, SFD use the correlation 



between intrinsic B — V color and the Mg2 line index in elliptical galaxies to convert from 
dust column density to reddening. It can be seen in Figure 6 of SFD that 90% of the 389 
elliptical galaxies they use for the B — V vs. Mg2 regression have E{B — V)sfd values of less 
then 0.1 mag. Moreover it is clear from their Figure 6 that the fit is excellent for values of 
E{B — V)sFD less than 0.15 mag, but for E{B — V)sfd > 0.15 mag the fit starts to diverge 
for the few points that exist. Galaxies with E{B — V)sfd > 0.15 mag {Ay > 0.5 mag) 
seem to follow a trend leading to reddening values in regions with high dust opacity being 
overestimated. SFD state that the slight trend in the residual is not statistically significant, 
but that may be due to the fact that there are so few points with E{B — V)sfd > 0.15 mag. 

Figure 6 in SFD shows only two points (galaxies) that have E{B — V)sfd > 0.3 mag. 
We use the point with E{B — V)sfd = 0.32 mag to asses the amount of overestimation by 
the SFD method in high reddening regions. Using Figure 6 in SFD and Ay = RyE{B — V), 
with Ry = 3.1, for an extinction of around 1 mag, the SFD method would overestimate 
the extinction, at 1.31 magnitudes. In fact, by comparing Ay^^^ to Ay^^^^, we find the 
overestimation by SFD can be even more than that: when Ay^^^^ is 1, Ay^^^^ is typically 
1.5. Thus we are convinced that SFD overestimates the reddening value to lines of sights 
where the extinction is more than 0.5 mag. Such overestimation is due to the fact that 
in the sample of 389 galaxies most of them have E{B — V) < 0.1 and very few have 
E{B — V) > 0.15. This results in an accurate conversion between 100 /im emission and 
reddening for regions with very low extinction {E{B — V) < 0.1), but a less accurate 
conversion where E{B — V) > 0.15 {Ay > 0.5 mag). 

One could argue that the SFD method overestimates the extinction because we 
overestimate the value of Ry, needed to convert E{B ~ V) to Ay. We believe that is not 
so, since several independent studies indicate that the value of Ry in the Taurus dark cloud 
complex is around 3.1 and that it stays constant through the region for lines of sight in 
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which Ay < 3 mag (Kenyon et al. 1994; Vrba & Rydgren 1985). Furthermore, Ry = 3.1 
gives very good agreement between Ay^^^^ and extinction obtained using the color excess 
of background stars for which spectral types are known (AG). 

Structure. In Section 2, we noted that there are certain places where the traces of 
^VissA ^^^ ^VsFD differ not only in value but also in structure. 

The most dramatic "structural" differences correspond to places where the cuts 
in Figure 1 intercept IRAS point sources. SFD remove point sources from the ISSA 
images before they use them, while AG do not. Both AG and SFD assume a single dust 
temperature for each line of sight. In AG's determination of Ay^^^^, the high 60 ^m flux 
from IRAS point sources transforms into dust with high color temperature, which in turn 
causes an artificial dip in the extinction (see Wood et al. 1994 and references therein for 
more on this effect). We can see that SFD have done a good job at removing point sources 
by noticing that where Ay^^^^ has artificial dips, Ay^^^ traces the real structure (also 
traced by the other methods in AG). Thus Ay^^^, unlike Ay^^^^, is free of any unphysical 
extinction due to IRAS point sources. 

Less dramatic, but more curious, structural discrepancies between the two methods are 
seen near the extinction peaks in Figure 1. Although Ay^^^ is in general larger than Ay^^g^ 
by a factor of 1.3 to 1.5, it is important to realize that in some peaks Ay^^^^ is larger by a 
factor less than 1.3 and in others it is smaller than Ay^^^^. In addition the gaussian-looking 
extinction peaks traced by Ay^^^ are wider than those traced by Ay^^^^. Gaussian fits to 
the peaks associated with Tau Ml and L1506 for both traces in cut 1 show that both peaks 
traced by Ay^^^ have a FWHM wider than those of the peaks traced by Ay^^^^ by a factor 
of 1.3. This suggests that in regions of steep extinction gradients the reddening obtained 
from the SFD map has a lower resolution than the quoted 6.1'. 

We believe that this "smearing" in the SFD maps is caused by the difference in 
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resolution between the ISSA flux maps (6.1') and the COBE-hmited color temperature 
maps (1.4°) employed by SFD. In Figure 2 we plot the 100 fira emission (/loo) versus 
declination for a section of cut 1, where the rise in /loo associated with Tau Ml and L1506 
are clearly seen. The plot shows the traces of Jioo obtained by SFD (Iioosfd) ^^^ ^J ^^ 
(Iiooissa) averaged over the ~ 10' width (in R.A.) of the cut. It can clearly be seen that 
both traces are essentially equal, and that the peaks have the same width and height. Thus 
it is after both the dust color temperature and Jioo have been used together to obtain 
the extinction, that the SFD method loses spatial resolution. Regions of steep extinction 
gradients are likely to have temperature gradients, but if the spatial resolution in the dust 
color temperature is ~ 1.4° then it is unlikely that these gradients would be detected. If the 
temperature gradients are not accounted for when calculating reddening, the result is areas 
with lower effective resolution. Therefore, the fact that SFD use a dust color temperature 
map with a spatial resolution a factor of 14 larger than the Jioo map results in a reddening 
map which does not have a constant spatial resolution of 6.1'. Unfortunately the DIRBE 
instrument has a limiting resolution of ~ 1°, so it is not possible to do better with data 
taken with that instrument. 

5. Conclusion 

We test the COBE/IRAS all-sky reddening map of Schlegel, Finkbeiner & Davis 
(1998) using the extinction study of Arce & Goodman (1999). In their study, AG use four 
different techniques to study the extinction in a region of the Taurus dark cloud complex. 
All four techniques give very similar results in terms of the value of Ay and the structure 
in the extinction. Thus the results of AG can be considered a truthful representation of the 
extinction in this region, and can be used to test the reliability of the SFD reddening map. 

We compare the extinction obtained by AG using ISSA 60 and 100 /im images with 
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the SFD reddening map. Our results show that in general the SFD method gives extinction 
values a factor of 1.3 to 1.5 larger than the extinction obtained by AG. We conclude that 
SFD overestimates the reddening value to lines of sights where Ay > 0.5 mag. We expect 
this overestimation is caused by the fact that in the sample of 389 galaxies used to calculate 
a conversion from dust column density to E{B — V), 90% of the galaxies have low reddening 
values {E{B — V) < 0.1) and very few (~ 4%) have high reddening {E{B — V) > 0.15) 
values. The lack of galaxies in high reddening regions results in inaccuracy in the conversion 
between dust column and reddening for lines of sight with E{B — V) > 0.15. This bias 
should be studied in other regions of the sky in order to see how general this trend is. For 
now we advise that caution be taken when using the SFD all-sky reddening map to obtain 
reddening (and extinction) values for regions with E{B — V) > 0.15 {Ay > 0.5 mag), as the 
value given by SFD could be larger than the real value. 

The behavior of Ay^p^ near steep gradients in extinction seems to be different from the 
overall "smoother" extinction regions. The difference between Ay^^^ and Ay^^^^ decreases 
near extinction peaks and in some peaks Ay^^^, even becomes smaller than Ay^^^^. In 
addition the peaks in Ay^^,^ are wider than the Ay^^^^ peaks. We attribute this behavior 
to the fact that SFD use a dust temperature map with a spatial resolution of a factor of 
14 larger than the 100 ^m intensity map used to obtain the reddening map. The poor 
resolution in the temperature maps results in a reddening map which shows less spatial 
resolution near extinction peaks. 

We thank David Schlegel, Douglas Finkbeiner and Mark David for their very useful 
comments on this work, and Krzysztof Stanek for pointing out the SFD paper to us, and 
for sharing his work on this topic. We also thank the referee for prompt and very helpful 
comments. 
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Fig. 1. — Plots of extinction vs. declination for cut 1 and cut 2 of AG. The first cut is 
centered on R.A. 4hl9m33s (B1950), and the second cut is centered on R.A. 4hl8m24s 
(B1950). The light solid line is the extinction obtained from the SFD reddening map. The 
dark solid line is the extinction obtained using the ISSA 60 and 100 /xm images, the light 
dotted line is the extinction from star counts and the dark dotted line is the average color 
excess technique from AG. The average 1 o random error for each of AG's extinction method 
is shown. The peaks in extinction associated with named dark clouds and IRAS cores are 
identified. Dips in the ^V/ssa trace, due to artificial drops in the extinction caused by the 
presence of IRAS point source in the images, are also shown. 
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Fig. 2. — A plot of 100 /iin dust emission versus declination for a part of cut 1. The light 
line is the emission obtained from the SFD 100 /im map. The dark line is the emission from 
the AG 100 yum map. The rise in 100 /xm emission associated with Tau Ml and L1506 are 
identified. 



